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Abstract: In A^ = 4 super- Yang-Mills theory at large N, large A, and finite temperature, 
the value of the Wilson-Maldacena loop wrapping the Euclidean time circle (the Polyakov- 
Maldacena loop, or PML) is computed by the area of a certain minimal surface in the dual 
supergravity background. This prescription can be used to calculate the free energy as a 
function of the PML (averaged over the spatial coordinates), by introducing into the bulk 
action a Lagrange multiplier term that fixes the (average) area of the appropriate minimal 
surface. This term, which can also be viewed as a chemical potential for the PML, contributes 
to the bulk stress tensor like a string stretching from the horizon to the boundary (smeared 
over the angular directions). We find the corresponding "hedgehog" black hole solutions 
numerically, within an S'0(6)-preserving ansatz, and derive part of the free energy diagram 
for the PML. As a warm-up problem, we also find exact solutions for hedgehog black holes 
in pure gravity, and derive the free energy and phase diagrams for that system. 
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1. Introduction 

There has recently been a lot of interest in the large- N thermodynamics of four-dimensional 
SU(N) gauge theories compactified on S s . Although they live on a space with finite volume, 
in the 't Hooft limit such theories have an infinite number of degrees of freedom and can 
therefore exhibit sharp phase transitions and spontaneous symmetry breaking. Much of the 
recent interest stems from the observation by Sundborg [[IJ and Aharony et al. Q that, even 
in the free limit, such theories can undergo a first-order phase transition as a function of 
temperature, and this transition shares enough features of the usual deconfinement transition 
in flat space to deserve the same name. Specifically, the free energy goes from being of order 
N° to A 2 , indicating liberation of the elementary colored fields, and the Polyakov loop (the 
Wilson loop wrapping the periodic imaginary time direction) becomes non-zero, indicating 
breaking of the Z n center symmetry. The strategy employed in that work was to compute the 
effective potential for a single mode of the system, namely the constant mode of the Polyakov 
loop over the three-sphere, by integrating out all other modes. From this effective potential, 
or off-shell free energy, one can read off the saddle points of the theory — stable, unstable, 
and meta-stable — and how they depend on temperature. Much of the subsequent work on 
weakly-coupled large N gauge theories on S 3 has continued to focus on the Polyakov loop 
and its effective potential g @, |, §, |, 



At large 't Hooft coupling, the easiest theory in this class to study is the conformal M = 4 
super- Yang-Mills theory, which is dual to type IIB supergravity on asymptotically AdSs x 
manifolds. At finite temperature, gravitational systems with a negative cosmological constant 
and asymptotically AdS boundary conditions are subject to a Hawking-Page phase transition 
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[jD| . This is a first-order transition, below which the dominant thermal state is AdS and 
above which it is a large black hole. (We briefly review this transition in Section 2.) From the 
gauge theory point of view, this is again a deconfinement transition []l4j], with the same two 
properties mentioned above: the free energy goes from being of of order TV to TV 2 , and the 
Ztv center symmetry breaks. The former effect is due to the fact that the black hole solution 
has a non-zero classical action, and this action carries an overall factor of 1/Gn ~ N 2 . The 
latter effect can be seen as follows. While the usual Wilson loop is not straightforward to 



compute using AdS/CFT (but see [15] for a proposal), one can compute a locally BPS version 
that includes a coupling to the six adjoint scalars (p 1 : 

W{C) = -J-TrM(C), M(C) = Pexp(f dsiiv^A^ + lvlO 1 ^ 1 ). (1.1) 
A Jc 

(Here = dx^/ds along the curve C, and 6 1 is an arbitrary unit six- vector.) This is the 
so-called Wilson-Maldacena loop, and it is given as e _Sdominant , where ^dominant is the saddle- 



point value of the action of a string world-sheet anchored on the curve C [16, 17]. In the 
absence of a -B-field, Sdominant is simply the area in string units of the minimal surface with 
those boundary conditions. (More precisely, this area is divergent due to the region near the 
boundary — a UV divergence from the boundary point of view — and must be renormalized.) 
Just like the Polyakov loop, the Wilson-Maldacena loop running around the imaginary time 
circle (the so-called Polyakov-Maldacena loop, or PML) is charged under the Z^r center. In 
the Euclidean version of AdS the imaginary time circle is non-contractible, so there actually 
aren't any world-sheets with the appropriate boundary conditions, implying that the PML 
vanishes and the center is unbroken. On the other hand, due to the existence of a horizon the 
imaginary time circle for a Euclidean black hole is contractible, leading to a finite value for 
Sdominant and hence the PML. The Lorentzian version of these statements is that the string 
ending on the horizon represents screening of the quark charge in the field theory; in the 
absence of a horizon the string has nowhere to end, and therefore inserting a single quark is 
forbidden. 1 

Much of the subsequent work on the thermodynamics of this system has focused on either 
the PML or the ordinary Polyakov loop, and their respective effective potentials [18, 1^, pOl 21, 



^2] , 23, 24, p5[ . In the papers [21, |24|, |25| , Wadia and collaborators studied effective potentials 
for the ordinary Polyakov loop based on phenomenological unitary matrix models. On the 
other hand, conjectural sketches of the free energy diagram for the PML, built around known 



properties of the saddle points, have appeared in various places, including [jj], [L9], |2C| |22|, 2q| . 



1 The discussion in this paragraph is strictly speaking valid only at infinite TV. At finite TV the system has a 
finite number of degrees of freedom and therefore no spontaneous symmetry breaking or phase transitions are 
possible. Correspondingly, the PML vanishes identically due to the Gauss law constraint. On the AdS side 
this is enforced by the integral over a certain 73-field Wilson line [Q. In order to see the effects described in 
this paragraph, therefore, it is necessary to take the large TV limit in the presence of a fixed but very small 
explicit breaking of the center. See Q for a detailed discussion of this point. In Section 3, we will define our 
order parameter in such a way that this issue does not arise, by also adding an anti-quark at the antipodal 
point on the S 3 . 
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Our purpose in this paper is to compute quantitatively the effective potential for the PML, 
using the gravitational dual theory. As we will now briefly explain, this will lead us to a 
problem in General Relativity that is of some interest in its own right. Due to the complexity 
of the gravitational theory, we will in this first pass only be partially successful, and part of 
the free energy diagram will remain a conjectural sketch. Based on our results, we will also 
sketch a phase diagram for the system in the presence of a chemical potential for the PML 
(roughly speaking a density of external quarks). 

The basic method we will use to compute the free energy diagram for the PML is due 
to York [^j; let us briefly review his work. He was studying gravity in a cavity at finite 
temperature (the thermodynamics of this system is closely analogous to that of asymptotically 
AdS spaces; the calculation described below was done in the asymptotically AdS case in 

|l9| , f20|). He considered sub-ensembles of the canonical ensemble in which the black hole 
horizon area was fixed. To do this, he added a Lagrange multiplier term to the Euclidean 
Einstein-Hilbert action. The saddle-point approximation to the free energy is then given by 
the action evaluated on the solution to the new equations of motion. The effect of the Lagrange 
multiplier term is to add to the stress tensor in the Einstein equation a term corresponding 
to a membrane wrapping the horizon. This stress-energy creates a conical deficit angle in the 
orthogonal (i.e. imaginary time-radial) plane. The result is a free energy diagram that clearly 
shows the main features of the system's thermodynamics: the three saddle points, namely hot 
flat space and large and small black holes, the thermodynamic instability of the small black 
hole, and the Hawking-Page transition (see figure 11 of ]28|1 ). Interestingly, the free energy 
turns out to be continuous even across the topology-changing transition separating the flat 
space and black hole regimes. 

The idea of applying the same method to computing the free energy diagram for the PML 
was discussed by Wiseman and the present author in the paper [28]. As above we should add 
a Lagrange multiplier term to the bulk gravity action, but this time fixing the area not of the 
horizon but of the orthogonal disc, spanned by the radial and imaginary time directions. Its 
effect is to add to the stress tensor a term corresponding to a fundamental string stretching 
from the horizon to the boundary. More precisely, to define the PML one must fix a point 
on the boundary 5 3 as well one as on the asymptotic 5 5 . To have the lightest possible 
mode as our order parameter, we take its (logarithmic) average over both spheres. (A precise 
definition of our order parameter is given in Section 3. It is morally similar to the one used by 
Aharony et al. [§, but here the averaging is done in a manifestly gauge-invariant, rather 
than gauge- fixed, manner.) We thus avoid explicitly breaking the system's SO (A) x 50(6) 
symmetry. Assuming that this symmetry is not spontaneously broken (which it is in some 
cases, as we will discuss below), the result is to smear the string source uniformly over both 
sets of angular directions. 

We are thus left with the problem of finding the solution to the Einstein equation for 
a black hole with a relativistic string stretching from the horizon to infinity and smeared 
over the angular directions. That problem was solved for pure gravity in four dimensions by 
Guendelman and Rabinowitz [29], who termed the solutions "hedgehog black holes". As a 
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warm-up problem, in Section 4 we generalize their solutions to arbitrary dimensions. The 
solutions are beautifully simple, as are the free energy and phase diagrams derived from them 
(see figures |2] and |3|). As in York's case, the free energy is continuous even as the topology 
changes from the black hole to AdS. Perhaps the most interest feature of the phase diagram 
is that the Hawking-Page transition curve ends at a critical point. 

The actual system of interest, supergravity on AdSs x S 5 , is discussed in Section 5. 
Thankfully, the symmetries of the problem remove most of the new fields; the only ones re- 
maining are the dilaton and a scalar representing the radius of the five-sphere. Thus, after 
gauge fixing, there are four degrees of freedom (the dilaton and the radii of the S , S 13 , and 
S 5 ) depending on one variable. Nonetheless the equations of motion are sufficiently compli- 
cated that they defy exact solution, and we had to resort to solving them numerically. The 
resulting free energy diagram is shown in figure ||. One unexpected feature is the existence 
of a (temperature-dependent) positive lower bound on the logarithm of the PML. However, 
it should be recalled that we have worked within an ansatz that preserves the SO (6) sym- 
metry of the boundary conditions. It is known that, at sufficiently high temperature, the 
conventional small black hole with that symmetry is unstable to perturbations that break 
it [p0| . It is also known that there exists at least one solution with lower free energy in 



which SO(6) is spontaneously broken, namely the ten-dimensional black hole [31|. This con- 
ventional solution should extend into a branch of ten-dimensional hedgehog solutions, which 
we conjecture continues all the way to unit PML. Finding these solutions, however, would 
present a challenging numerical problem, and no attempt is made to do so here (even the 
conventional ten-dimensional black hole is not known analytically or numerically, outside of 
the limit that its horizon is much smaller than the AdS radius In the phase diagram 

for the system with a chemical potential for the PML, the Hawking-Page transition should 
extend to a first-order transition separating five- and ten-dimensional black hole phases. We 
sketch a possible such phase diagram in figure |8|. 

For completeness we have also done the analogous calculations for the gauge theory on 
R 3 , which can be considered as the high-temperature limit of the theory on S 3 (subsection 
5.2). 

Let us close this introduction with two general comments. First, one might ask what hap- 
pens if one retains the entire matrix M ( |l.l| ) whose trace is the PML, rather than integrating 
out everything but the PML itself. In many instances such matrix models capture the inter- 
esting physics of a system more directly than is revealed by computing the effective potential 
for the trace alone. Finding the matrix model governing M is equivalent to computing the 
effective action for the collection of all possible traces, 

W ninana ... = ^ Tr(M ni Mt^M^ • • • ) , (1.2) 

where (ni,ri2,n3, . . . ) is an arbitrary collection of (positive or negative) integers (recall that 
M is a general matrix, neither Hermitian nor unitary). Unfortunately, except for the simple 
traces of the form W n = Tr M n /N and Tr M' n , it is not clear how to compute such observables 
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using AdS/CFT. Furthermore, even those simple traces are not independent observables in the 
supergravity limit, since they obey the identity W n = (up to a' and 1/N corrections), 
since the minimal surface for the n-fold cover of a curve is simply the n-fold cover of the 
basic minimal surface (putting aside the special case where the fundamental group of the 
bulk topology has a finite part). Thus it is not possible, by adjusting the supergravity 
configuration, to independently fix W n for different n. To summarize, in the supergravity 
limit the effective potential for W is "all there is" . 

A second comment is that, in the papers [18, 19, [H], 21, 22, 24], one of the motivations 
for going off-shell was to study the Hagedorn transition of string theory, which on-shell is 
hidden behind the Hawking-Page transition. In particular these papers conjecture that the 
Hagedorn transition occurs when the thermal AdS and small black hole (or possibly ten- 
dimensional black hole) saddle points merge. Unfortunately, since we work within the gravity 
approximation, our free energy diagram is too crude to address this stringy issue. However, 
as a matter of principle our general method for going off-shell should be applicable to the full 
string theory; it would be interesting to see whether it can be applied in practice. 



2. Review: Hawking-Page transition 

The purpose of this section is to remind the reader of the basic facts about the Hawking-Page 
transition 13], while setting up the framework in which we work. 



We consider four-dimensional SU (N) Af = 4 super- Yang- Mills theory, compactified on a 
unit three-sphere. In the regime of parameters 1 A ^ iV 2 , where A = gy^N is the 't Hooft 
coupling, the theory is best described by type IIB supergravity (UV completed by type IIB 
string theory) on asymptotically AdSs x 5 5 manifolds. The curvature is supported by N units 
of five- form flux on the S 5 . Stringy effects are suppressed by 1/A and quantum effects by 
1/N; more precisely, setting the S 5 and AdSs radii to one, the ten-dimensional gravitational 
coupling constant is Gn = ^7r 4 A r ~ 2 while the string length is a' 1 / 2 = A -1 / 4 . 

The partition function for the canonical ensemble at temperature T is defined by the 
Euclidean path integral on S 1 x S 3 , where the imaginary time circle has circumference 1/T 
and anti-periodic boundary conditions for the fermions. It is computed by the gravitational 
path integral over Euclidean manifolds that are asymptotically AdS^ x S 5 , where AdS^ is 
Euclidean AdSs with "imaginary time" coordinate periodically identified so that its conformal 
boundary is the same S 1 x S 3 on which the field theory lives. Again, at large N the saddle 
point approximation is valid, and we have (to lowest order in 1/N) 



e 



-F/T _ ~ g — ^dominant (2 \) 



where Sdominant i s the minimum value of the Euclidean action among the saddle points. 
Actually, given the boundary conditions, on any solution this action will be divergent due 
to the infinite volume near the boundary (an ultraviolet divergence from the field theory 
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viewpoint), and must therefore be regularized 2 and made finite by a divergent subtraction; 
conventionally the counterterm is chosen so that AdS^ x S 5 has vanishing renormalized action. 
Since the action comes multiplied by an overall 1/Gn ~ N 2 , the free energy is of order TV 2 , 
unless the action of the dominant saddle point vanishes. In the latter case the one-loop 
contribution, which is of order N°, is the leading term; it represents the contribution from 
the thermal gas of gravitons. 

At any temperature, AdSg x S 5 is a saddle point, and at temperatures below T cv n — 
8 1 / 2 /(27r) it is the only one. Above T cr ; t there are at least two others, namely the small and 
large AdS-Schwarzschild black holes. Both solutions take the form of a product of S 5 with 
the following asymptotically AdS^ Einstein metric: 



ds z = f(r)dr z + f{r)~"dr z + r z dn z 3 , (2.2) 

/(r)=r 2 + l-4, t~t+-, r h <r<oo, (2.3) 
r z 1 

where the horizon radius is the largest zero of f(r), and is fixed by the requirement 
of smoothness at the horizon, /' (i"h) = 47rT. None of the supergravity fields aside from 
the metric and five-form are excited. Eliminating ji from the two equations /(r^) = and 
f'( r h) = 47rT, we arrive at a quadratic equation for r^. For T > T crit , this equation admits 
two solutions, one with < 2 -1 / 2 , the other with > 2~ 1//2 ; these are the small and large 
black holes respectively. Their actions are 

N 2 

S = — r 2 (l - r 2 ) . (2.4) 

Any black hole with > 1 is therefore thermodynamically dominant over thermal AdS; 
this is never the case for the small black hole, but is the case for the large black hole when 
r>r H P = 3/(2vr). 

Besides never being thermodynamically dominant, the small black hole has a negative 
eigenmode of the Lichnerowicz operator |32| (similar to the Schwarzschild solution's Gross- 
Perry- Yaffe mode |33|). This is related to the fact that its Lorentzian version represents a 
black hole in an unstable thermal equilibrium with the heat bath, due to its negative specific 
heat. Furthermore, for less than around 0.426, which translates to temperatures above 
around 3.20/(27r) (slightly higher than Thp), the small black hole suffers from a second, 
dynamical instability, namely a Gregory-Laflamme instability on the S 5 [ 30 1 . 

The large and small AdS-Schwarzschild black holes are uniform on the 5 5 , and therefore 
preserve the boundary conditions' SO (6) symmetry. At sufficiently high temperatures, there 
are other saddle points that spontaneously break that symmetry. These presumably include 
a black hole that is spread out non-uniformly over the S 5 , and a ten-dimensional black hole 
[ j3ll l that is localized on the S 5 and breaks the SO (6) down to SO (5); its horizon topology 



2 The simplest regularization procedure is to impose a finite Dirichlet boundary condition on the metric. We 
require the induced metric to be Sri^t) x Sr x &i> where the subscripts indicate the radii of the respective 
spheres. R is taken to infinity to remove the regulator. 
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is S 8 , rather than S 3 x in the case of the AdS-Schwarzschild black holes. At high tem- 
peratures, where its horizon is very small compared to the AdS radius, this is essentially a 
ten-dimensional Schwarzschild black hole, and is therefore presumably dynamically stable but 
thermodynamically unstable. Exact solutions for these saddle points (or others with more 
exotic horizon topologies) are still unknown, as are their precise properties, such as at what 
temperatures they appear. 



3. Order parameter and off-shell free energy 



The Hawking-Page phase transition in AdSs x S 5 described in the previous section is, from 
the point of view of the boundary gauge theory on S 3 , a deconfinement transition [ 14 1 . There 



are two features which qualify it for this name. First, as explained above, the free energy goes 
from being of order TV to iV 2 . Second, the Z^v center symmetry breaks. In this section we 
will study this breaking by defining a certain order parameter for the symmetry. We will then 
explain how, in the rest of the paper, we will compute the off-shell free energy as a function 
of that order parameter. 

We start by defining the Polyakov-Maldacena loop (PML). This is a function on S 3 x S 5 . 
Using r and 6 to denote the coordinates on the S 1 and S 3 respectively, and 9 1 to denote a 



unit vector in R , we have [16, 17 1 



W(9,9) = -^TiPexp jdr (iA T {r,9) + #VM)) . (3.1) 

Here A is the gauge field and tfi are the six adjoint scalars. W has unit charge under the Zjv 
center symmetry; it could therefore serve as an order parameter, except that, for the very 
same reason, its expectation value vanishes identically (no spontaneous symmetry breaking 
can occur in a system with a finite number of degrees of freedom 3 ) |ji~4"|| . Therefore we multiply 
W by the PML for an antiquark, which has charge —1 under the Zjv, that we choose to 
insert at the same point on the S 5 but at the antipodal point —9 on the S 3 : 

w(9,9) = W(9,9)W j (-9,9). (3.2) 

Although neutral under the Zat center of the gauge group, w nonetheless serves as an or- 
der parameter; its expectation value is of order iV -2 when the symmetry is unbroken, and 
otherwise of order iV . 4 We have 

w = e~ f ™ /T , (3.3) 



3 An equivalent way to see that its expectation value vanishes is to use the Gauss law constraint; the charge 
of a single external quark cannot be neutralized by the elementary quanta of the theory, which all transform 
in the adjoint representation, to produce a singlet under the global part of the gauge group. 

4 This fact is explained in detail in |B| . It can be understood roughly from the fact that large N factorization, 
which would have implied (w) = (W)(W^) + 0(N~ 2 ), fails only because in the path integral we must sum 
over gauge configurations that are related by the action of the center symmetry; this sum sets (W) to zero but 
not (w). Aside from this subtlety, the operators we deal with become classical in the large N limit, so we will 
generally simply refer to their "values" rather than their "expectation values" . 
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where f„q is the free energy cost of adding to the system an external quark and anti-quark 
(charged under the scalars as well as the gauge field). Because of the self-energies of the quark 
and anti-quark, f qq - suffers from an ultraviolet divergence and must be renormalized. 

In the gravity dual, the external quark and anti-quark source fundamental strings that 
extend into the bulk (with opposite orientations, i.e. if we put an arrow on the string, then 
the arrow would point away from the quark and towards the anti-quark). If we are at large 
N, then as in the previous section we appeal to the saddle-point approximation and imagine 
that the bulk is described by a fixed classical supergravity solution. If we are also at large 
A, then the quantum fluctuations of the world-sheet are also suppressed, and to a good 
approximation the string world-sheet lies on a fixed surface in the bulk geometry, namely 
the solution to its classical equations of motion (if there are several solutions, then the one 
with the smallest Euclidean action). The quark and anti-quark world-lines on the boundary 
provide the boundary conditions for this surface, along with the value of 8, which fixes the 



asymptotic position of the string on the S . We have 17 1: 



w = j^e Sdominant , (3.4) 

where we use s for the Euclidean world-sheet action. (The factor of N~ 2 arises from the same 
factor in the definition of w.) s dominant is divergent due to the infinite area of the world-sheet 
near the boundary (reflecting the ultraviolet divergence of f qq ), and must be renormalized, 
similarly to the bulk action S. 

The world-sheet action s is, for our purposes, made up of two terms. The topological 
(Fradkin-Tseytlin) term contributes a factor of N 2 to e~ s if the world-sheet is topologically 
two discs, and N° if it is a cylinder (just as in the usual 't Hooft genus counting). Therefore 
a saddle point with two discs, if it exists, will always dominate in 't Hooft limit over a 
cylinder. The Nambu-Goto term, which is the (string-frame) area of the world-sheet in string 
units, requires the world-sheet to be the minimal surface within each topological sector. (If 
the supergravity background has a -B-field, then it also needs to be taken into account in 
minimizing the world-sheet action. Non-zero Ramond-Ramond fields — which are of course 
present in the case at hand — couple only to the world-sheet fermions and therefore should 
not affect the classical solutions.) 

Let us first consider the case where the bulk metric is the Euclidean black hole solution 
( |2.2D . The dominant world-sheet configuration has two connected components, which wrap 
the disc and are constant on the S s at the points 9 and —6 respectively, and also constant on 
the S 5 at the point 6. The area of the disc is easily read off from the metric fl2.2|) [B2]: 



A ct + ± r d r = -^, (3.5) 

where we have added an ^-independent counterterm {—R/T in the regularization described 
in footnote |2|) . We thus find 

f ™ = -^ = -—- (3 - 6) 
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On the other hand, when the bulk metric is thermal AdS, the Euclidean time-circle is 
non-contractible, so there are no world-sheet configurations that are topologically two discs. 
The dominant one is instead a cylinder, with the string connecting the quark to the anti-quark 
through the center of AdS, and wrapping the S 1 (again at a fixed position on the S 5 ). This 
world-sheet has vanishing renormalized area, so we have 

w = ■ ( 3 - 7 ) 

In the rest of the paper, we will be interested in fixing a particular value for the PML, 
and integrating out all other degrees of freedom. More specifically, we will take as our order 
parameter the average value of f q q over the three- and five-spheres (times — 7rA -1 / 2 ): 5 

The reason for taking the average over the spheres is to retain the lightest possible mode 
when we integrate out the other modes. In this way, we avoid explicitly breaking the theory's 
SO(4) x SO(6) global symmetry. The factor of A -1 / 2 is for convenience, in order to have 
a quantity which is finite in the limit A — > oo. We have, for example, $ = for AdS- 
Schwarzschild and $ = — 2nT\~ 1 / 2 In N for pure AdS (the latter value goes to — oo in the 
strict 't Hooft limit; on the other hand it goes to zero if we take N to infinity holding, say, 
<7ym fixed, which is also well within the regime of validity of supergravity) . Those solutions 
both have SO (4) x SO (6) isometries which dictate the position of the classical string world- 
sheet, and therefore <J> reduces to a local functional of the metric; we will lean heavily on this 
convenient property in our computations below. In general, without those isometries, $ is a 
non-local functional of the metric since one must solve the string equations of motion. 
We now define the effective potential (or off-shell free energy) F(<&) by 



e 



F(*o)/T = / _ $ ) e -m . (3.10) 



Here <fi includes all the degrees of freedom of the system — in our case, all modes of all the 
supergravity fields. We then have 



e 



W = z= d®e- F W' T . (3.11) 



5 It is interesting to compare our order parameter with the one employed by Aharony et al. [y] for their 
weak-coupling calculations: 



1 Tr P exp j dr i J »^- flx ed (r , Q) , (3.8) 

where they have picked a particular gauge. Gauge-fixing is necessary for them in order to make the average 
of A T over the three-sphere a sensible (i.e. gauge-invariant) observable. Thus one difference between our order 
parameter and theirs is that we take this average after computing IniV -1 TrPexp^dr, whereas they take it 
before. Another difference is that they do not include a coupling to the scalars (therefore they have no need 
to average over S 5 ). 
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In the thermodynamic limit (in our case, the large N limit) this last integral is evaluated using 
the saddle point method, i.e. simply by finding the global minimum of F(<&), which represents 
the thermal equilibrium; other local minima represent metastable equilibrium configurations. 

Our interest will be in evaluating the right-hand side of ( |3.10| ), which we also will do by 
a saddle-point method. In order to restrict the integral to configurations with the given value 
of $>[</>], we add to the supergravity action a Lagrange multiplier term, 

SsugraM -> SsugraM - £ m<P] - *o) • (3-12) 

We then find the solutions to this modified action, and evaluate the action on them. The 
effect of the new term on the equations of motion will be discussed in detail in the next two 
sections. 



4. Warm-up: pure gravity 

In this section we will solve the problem explained in the last section for a simpler theory 
than type IIB supergravity, namely pure gravity with negative cosmological constant. We 
will see that this system is simple enough to allow an exact solution. 

It turns out to be just as easy to work in a general spacetime dimension D = p + 2 
(p > 1). As above we set the AdS radius to one. Again, we work in the canonical ensemble 
at a temperature T, so we consider Euclidean metrics that are asymptotically AdSp +2 > that 
is, whose conformal boundary is S\,, 2n T) x ^T- These boundary conditions have an SO (2) x 
SO(p + 1) symmetry, which we assume is shared by the metric; in other words we look for 
static, spherically symmetric solutions. We thus take the following metric ansatz: 

ds 2 = g TT (r)dT 2 + g rr (r)dr 2 + gnn(r)dnl , r~r+^, lim = 1 . (4.1) 

T r^oo g nQ 



In this context, our order parameter ( p.9p becomes 

/dP 6 f°° 
— drdr (g rr g T r) 1/2 = *«* - / dr {g rr g TT ) 1 ' 2 . (4.2) 

Varying the Lagrange multiplier term — — §q)/T with respect to the metric gives a 

contribution to the stress tensor which is that of a relativistic string of tension k extended 
in the r and r directions and smeared over the angular directions. The Einstein equation is 
now: 

G\ = G\ = \p(p + 1) - , G% = \p(p + 1) . (4.3) 

The first term on the two right-hand sides is the cosmological constant. On the right-hand 
side of G T T and G r r , the string tension k is multiplied by the gravitational coupling and 
divided by the proper area of the p-sphere (since the string tension is spread out over that 
area) . 
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It is well known (see |p4| for a discussion) that, in the presence of spherical symmetry, 
if the stress tensor obeys T T T = T r r and T T r = T r T = 0, then the solution to the Einstein 
equation can be written in the Schwarzschild form, 

ds 2 = f(r)dr 2 + /(r)- 1 * 2 + r 2 d^ 2 p . (4.4) 

The function / is easily found, and is a remarkably simple generalization of the AdS-Schwarz- 
schild solution: 

f(r) = r +1 — x T . (4.5) 

As in the case of AdS-Schwarzschild (Section 2 above), \x is fixed by the requirement of 
smoothness of the metric at the horizon, which requires f'( r h) = 47rT, where the horizon ra- 
dius is the largest root of /. The solution ( |4,4[ ), ( [4.5| ) was found (for p = 2) by Guendelman 



and Rabinowitz [29], who termed it a "hedgehog black hole". 

At long distances, where the gravitational field is weak, the hedgehog black hole metric 
can be derived from Newtonian reasoning as follows. As usual in the weak field regime we 
have f(r) = 1 + 2V(r), where V is the Newtonian gravitational potential. Thanks to the 
spherical symmetry, the latter is given by 

where M(r) is the total mass contained within a ball of radius r. In the case of AdS- 
Schwarzschild, M(r) includes a negative contribution proportional to r p+1 from the vacuum 
energy, plus an arbitrary constant representing a mass located at the origin. In the present 
case there is also a contribution nr due to the string. The fact that this Newtonian reasoning 
gives the correct exact solution to the full Einstein equation can be traced to the facts that, 



for a metric of the form (4.5), the Einstein tensor (with one raised index) is linear in /, and 
the stress tensor (with one raised index) is independent of /. 6 

We wish to compute, for the family of metrics obtained above, the relation between the 
order parameter <I> and the free energy F, which is T times the (renormalized) Euclidean 
action. $ is rather trivial to compute, and in fact doesn't depend at all the particular form 
of /; as in the case of AdS-Schwarzschild, we have 

<D = r h . (4.7) 

The free energy can be obtained by directly computing the Euclidean action in an appro- 
priate regularization and subtraction scheme (such as the one discussed in footnote |2|) . The 
calculation is straightforward but rather tedious. A simpler but less direct approach, which 
yields the same result, is to first find the relation between $ and k, and then use the relation 

s- 

6 The Newtonian reasoning continues to give the correct solution if we also include a spherically symmetric 



electric field; we simply add the term Gn? 2 /r 2p 2 to the right-hand side of (4.5), where q is the electric charge 
of the hedgehog black hole. 
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to find F as a function of <I>. The relation between <3? and k can be found by combining the 
equations /(r^) = and /'(rh) = 47rT to eliminate fi, and then using ( |4.7| ); the result is: 



167tGnk 



(p + l)$ p - 47rT$ p - 1 + (p - l)$ p " 



from which we obtain 



(4.9) 



(4.10) 




The integration constant can be fixed, for p > 2, by noting that the solutions go over in 
the limit $ — > to AdS p +2, which by convention has vanishing free energy. For p < 2, it is 
necessary to resort to the direct calculation of the Euclidean action. 

In figure |l] we plot F versus for p = fg n 
1, and in figure ^ both k and F for p = 3, 
at various temperatures. All of the im- 
portant thermodynamic features are clear 
from these figures: for p = 3, the appear- 
ance of the large and small black holes at 
T cr it, the instability of the small black hole 
and stability of the large one, and the fact 
that the large black hole becomes domi- 
nant at Thp; for p = 1, the absence of 
a small black hole and the fact that the 
large black hole exists even at zero tem- 
perature (although it is not thermodynam- 
ically dominant there). 

Let us focus on the case p = 3, and 
ask where AdSs belongs on the free energy 
diagram. (Strictly speaking, since we are dealing here with a toy model, this is an artificial 
question; however, let us pretend for the moment that we are dealing with a real string theory.) 
As mentioned in the previous section, for AdS we have <3? = — 2nT\~ 1 l' 2 In N, whose value 
at large N depends on how one takes the limit. The conventional 't Hooft limit (N — > oo, 
A fixed) yields $ = — oo. The free energy diagram then seems to suffer from a horizontal 
discontinuity, with no solutions connecting AdSs at $ = — oo to the hedgehog black holes at 
$ > (the hedgehog black holes with $ < are nakedly singular). However, we should recall 
that in the 't Hooft limit the string length A -1 / 2 remains fixed. In particular, hedgehogs with 
$ of order A -1 / 2 will experience large stringy corrections near the horizon (recall that = <£). 
In other contexts, such corrections are known to smooth out topology-changing transitions; 
one may therefore expect that, once they are taken into account here, the hedgehogs will 
extend all the way to $ = — oo and connect smoothly to AdS. Less speculatively, one can 
choose to take the large limit in such a way that the stringy corrections disappear at the 
same time as the quantum corrections, for example by keeping gy M fixed. This puts AdS at 



Figure 1: F versus $ in three dimensions (p = 1), 
at 4irT = 0,1,2,3 (top to bottom). The minimum 
of each curve is a BTZ black hole (at T = the 
minimum is at $ = and represents the zero-mass 
black hole) . The third curve is at the Hawking-Page 
temperature. 
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Figure 2: k and F versus <£> in five dimensions (p = 3), at (2irT) 2 — 5, 6, 7, 8, 9, 10 (top to bottom). 
From the bottom, the second curve is at Tjjp, the third at T cr it, and the fifth at T^. it (where the curves' 
inflection points merge and disappear). 

$ = 0, and gives a continuous free energy diagram. This is the point of view we will adopt 
in the rest of the paper. 

If we consider k to be a chemical potential rather than a Lagrange multiplier, in other 
words if we change the action by —K<&[g]/T rather than —n(&[g] — $q)/T, then we can use 
the above free energy diagrams to derive a phase diagram in the k—T plane. The effect of 
the chemical potential is to add a density of external "quarks", which source strings that 
extend into the bulk and end on the horizon. Figure |3| shows the phase diagram for the five- 
dimensional (p = 3) case. To understand this diagram, let us work at fixed temperature and 
study the effect of varying k. Above the Hawking-Page temperature, the chemical potential 
has no qualitative effect, since the system is already in the large black hole state (it does 
make the horizon larger — the strings effectively pull on it). However, below the Hawking- 
Page temperature, where in the absence of a chemical potential there is no black hole, the 
strings force a horizon to open. This can be seen from the free energy diagram, which is 
quadratic in the neighborhood of the origin. Adding a chemical potential creates a local 
minimum at <3? ~ 4Gn k / (Sir) , which, for sufficiently small k, is the global mininum. There is 
thus a second-order transition at k = 0, from the thermal AdS state to a "tiny black hole" 
(not to be confused with the small black hole, which is thermodynamically unstable and 
therefore never dominant). For temperatures between Thp and TLu = 6 1//2 /(2-7r), there is a 
further, first-order transition between the tiny and the large black holes, at 




(4.11) 
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u 1 kG n 
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4-y6 

Figure 3: Phase diagram in the n—T plane for the five-dimensional system (p = 3), whose free energy 
diagram is shown in figure |[ The dashed line represents a second-order transition at an infinitesimal 
value of k separating thermal AdS from a tiny black hole. The solid curve represents a first-order 
transition separating the tiny and large black hole phases. This is an extension of the Hawking-Page 
transition and ends at a critical point,. 

As we lower the temperature past T^ rit , the tiny and large black hole saddle points merge, 
and the phase transition curve separating them ends at a critical point. 

The same calculation can easily be done in the Poincare patch of AdS, where the boundary 
"gauge theory" lives on R p rather than S p . The resulting hedgehog black brane solution is 

ds 2 = /(r)dr 2 + /(r)-W + rV p , f(r) = r 2 - - JL. . (4.12) 

The effective potential for <& is 

16vrG N F = p<S> p+1 - 4vrT$ p . (4.13) 

(Here k and F are the string tension and free energy per unit volume on the boundary.) The 
Hawking-Page transition in this system occurs at zero temperature and is second-order. 

5. Supergravity 

5.1 Gauge theory on S 3 

In this section, we extend the calculation of the previous section to type IIB supergravity 
on AdSs x S 5 . The boundary conditions for the problem have an SO (2) x 50(4) x 50(6) 
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symmetry. For the time being, we will work within an ansatz for the supergravity fields that 
respects this symmetry; below, we will consider the issue of spontaneous symmetry breaking 
(particularly of the SO (6) factor). The only supergravity fields allowed by these symmetries 
are the following: the metric, which must be of the form 

ds 2 = g TT {r)dT 2 + g rr {r)dr 2 + e 2a{r) d^l + e 2a(r) dnj , (5.1) 

(this is the string-frame metric); the dilaton <J?(r) (not to be confused with our order param- 
eter); the Ramond-Ramond scalar, which is not sourced by the Lagrange multiplier term and 
can consistently be set to zero; and the five-form field strength, which, taking into account 
the flux quantization on the S 5 , must be of the form 

F 5 = 16irNa' 2 (d 5 9 + *d 5 #) , (5.2) 

where d 5 9 is the volume form on the unit S 5 . The five- form thus does not carry any inde- 
pendent dynamical degrees of freedom, and we have in total five degrees of freedom minus 
one gauge degree of freedom. Within this ansatz, the supergravity equations of motion can 
be derived from the following action: 

S= ~^T l dr (9rr9rr) 1/2 e- 2 ' t ' 

x (R(2) + 9 rr (4(9 r 0) 2 - 3(d r a) 2 - 5(d r a) 2 ) + 6e" 2<T + 20e" 2S - 8e 2 *+ 3CT - 5 ^ , (5.3) 

where Rm) is the Ricci scalar of the two-dimensional metric g TT dr 2 + g rr dr 2 , and = — 
|<7 — \a is the two-dimensional dilaton. (When directly computing the on-shell value of the 
Euclidean action, one must also take into account a counterterm and boundary term.) The 
order parameter is, as in ( [4,2| ), 

<%] = - j ' dr (g rr g TT ) l l 2 . (5.4) 

When added to the action, the Lagrange multiplier term — k(<£[<7] — §q)/T sources the metric, 
which sources eft, which in turn sources a and a. Therefore all the above fields are active in 
the solutions. 

After gauge fixing, the equations of motion derived from the action S — K{f&[g\ — $q)/T 
reduce to a system of four coupled ordinary differential equations. The boundary conditions 
on the horizon are determined by the requirement of smoothness of the metric and scalar fields. 
On the outer boundary, at r = oo, we have asymptotically AdS<^ x S 5 boundary conditions 
(with Dirichlet boundary conditions on the dilaton); in practice these boundary conditions 
were imposed by matching onto an asymptotic solution at a finite value of r. The resulting 
boundary value problem was solved numerically using a shooting method, implemented in 
Mathematica using a combination of the NDSolve and FindRoot routines. 

The results are shown in figure || in the form of plots of k versus at the same six 
temperatures as in figure |2|. The points where the curves intersect the horizontal axis represent 
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Figure 4: k versus $ for the same set of temperatures as shown in figure namely (from top 
to bottom) (2ttT) 2 = 5,6,7,8,9,10. 2ttT = 8 1/2 is the critical temperature, while 2ttT = 3 is the 
Hawking-Page temperature. (The curves for some higher temperatures are shown in figure ^[) The 
solutions corresponding to the points labelled a, 6, c, d are plotted in figure |^. 

conventional black holes. Above the critical temperature T cr jt = 8 1 ^ 2 /(2tt) (bottom two 
curves), the curves cross the axis twice, representing the large and small black holes. At the 
critical temperature (third from bottom) it is tangent to the axis, while below the critical 
temperature (top three curves) they do not intersect it at all. The major new feature compared 
to the same plot for pure gravity (left side of figure |2|) is that, whereas the latter curves go 
through the origin, the curves for the supergravity system seem to be "repelled" by the 
origin. At each temperature there is a positive minimum value of <3? below which there are 
no solutions, and above which there are two. Thus it would seem that we have failed in our 
attempt to fix the value of the order parameter, at least to certain values. We will discuss a 
possible resolution to this issue in the next subsection. It is also interesting to note that the 
upper branches of the curves quickly merge with each other. Inspection of the corresponding 
solutions shows that those with the same but different T on that branch are essentially 
identical outside of a neighborhood of the horizon, and that neighborhood becomes smaller 
and smaller as $ increases. A few examples of solutions are plotted in figure [|. 

The free energy can be computed from the curves of figure || by using the relation 
dF/d& = k. 7 It is plotted in figure |6[ For each temperature, we have only plotted the 
branch of solutions with the lower value of k, since these clearly also have the lower value 

7 The integration constants were fixed as follows. First, the large black hole at the Hawking-Page tempera- 
ture (2ttT = 3, $ = 1) has vanishing action, fixing the integration constant for that temperature. The other 
temperatures can then be fixed by matching in the region where the curves for the different temperatures 
merge, namely along the upper branch of solutions in the $-/c plane. 
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Figure 5: A selection of supergravity hedgehog black hole solutions, corresponding to the labelled 
points in figure |. a: (2vrT, $, n/N 2 , F/N 2 ) = (10 1 / 2 , 0.37, O.f 5, 0.04); 6: (fO 1 / 2 , 0.83, -0.30, -0.04); 
c: (f 1 / 2 , 1 .14, 0, — O.f 0) (this is the conventional large black hole at this temperature); d: 
(5 1 / 2 , 0.55, 0.33, 0.09). The horizontal axis is proper radial distance from the horizon. The curves 
marked S 1,3 ' 5 represent the proper radii of the respective spheres, and the one marked e* represents 
the local string coupling (relative to its asymptotic value g s ). In case c (the conventional black hole) 
both the dilaton and S 5 radius are constant. 

of F. The curves are qualitatively similar to the case of pure gravity (right side of figure 
H), and again exhibit the major thermodynamic features of the system: the appearance of 
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Figure 6: Free energy diagram for the same set of temperatures as in previous figures. Only the 
lower branch of solutions is plotted, since they obviously have lower free enegy at a given value of $ 
than the solutions on the upper branch. 



the large and small black holes at T cr jt, the thermodynamic stability of the large black hole 
and instability of the small one, and the fact that the large black hole becomes dominant at 
Thp- The major difference with the case of pure gravity is, again, that each curve ends at a 
minimum value of <£. 

Our discussion so far has neglected the possibility of spontaneous breaking of the sys- 
tem's SO(2) x SO(4) x SO(6) symmetry. Searching for solutions to the full ten-dimensional 
supergravity equations of motion, without imposing those symmetries, would present a chal- 
lenging numerical problem, and no attempt was made to do so here. (Aside from the increase 
in dimensionality, a complication that arises when the symmetries are not imposed is that, 
as mentioned in Section 3, &[g] becomes a non-local functional. To include the effect of the 
Lagrange multiplier, it is therefore necessary to solve simultaneously the string's equations 
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Figure 7: Conjectural sketch of the complete free energy diagram, at the representative temperature 
Tjjp, when the possibility of spontaneous 50(6) breaking is included. The upper curve is the same one 
shown in the previous figure, and represents 5'0(6)-preserving hedgehogs. The points marked a and 
b are the five-dimensional small and large black holes respectively. The lower curve is a conjectural 
sketch of a new branch of S'0(6)-breaking hedgehogs, which range from <f> = (thermal AdS) to 
$ = $*, where they connect continuously onto the upper curve. The maximum of this new branch (c) 
is a conventional black hole (since it has k = 0), and is presumably the ten-dimensional black hole, 
which is thermodynamically unstable but otherwise stable. The value of could be determined by a 
linearized perturbation analysis of the S'0(6)-preserving hedgehogs; the rest of the curve would require 
a difficult computation (a warm-up problem would be to find the conventional ten-dimensional black 
hole solution c). Somewhere between the point c and the new branch of hedgehogs would undergo a 
topology-changing transition, from ten-dimensional black holes localized on the S 5 to five-dimensional 
black holes that are non- uniform over the S 5 . 



of motion and the gravitational ones. Up to now we have only had to solve the gravitational 
equations, since the symmetries restrict the string world-sheet to lie in the r-r plane at fixed 
9 and 9.) However, some information can be obtained from the known properties of the 
conventional black holes in this system, which were reviewed at the end of Section 2. In 
particular, the large black hole is always locally stable, whereas the small black hole suffers 
from a Gregory-Laflamme instability on the S 5 at temperatures above Tgl ~ 3.20/(27r) |30t| . 
While we have not done the linearized perturbation analysis for our hedgehog solutions (it 
should not be too difficult in principle), the simplest picture consistent with the above facts 
is the following: at each temperature, there is a critical value $*(T) such that, for > ^(T) 
the hedgehog is stable, and for $ < 3>*(T) it is unstable against SO(6)-breaking perturba- 
tions; furthermore, for T in the range T crit < T < T GL , $*(T) < $ small black hole, while for 
T > T G l, $smaii black hole < $*{T) < $i a rge black hole- For $ < there should exist a branch 
of SO (6)-breaking hedgehogs with free energies less than the S'0(6)-preserving ones plotted in 
figure ||, which continuously connects onto those solutions at This new branch would in- 
clude both ten-dimensional black holes localized on the S 5 , and five-dimensional ones that are 
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non- uniform on the S 5 . The ten-dimensional conventional black hole should correspond to a 
local maximum on this new branch (maximum because that black hole is thermodynamically 
unstable). Finally, we conjecture that that branch extends all the way down to $ = 0, where 
it continuously connects onto thermal AdS, as occurred for pure gravity; the 50(6) symmety 
is restored as the tiny 10-dimensional black hole disappears. Figure |?] shows a sketch of this 
conjectured branch of 50(6)-breaking hedgehogs, superimposed on the SO(6)-preserving free 
energy diagram at the Hawking-Page temperature. 

2nT 



large five-dimensional black hole 




kG n 



Figure 8: Conjectural sketch of the phase diagram in the presence of a chemical potential k for the 
PML. 

If the free energy diagram sketched in the previous paragraph and in figure is correct, 
then we can derive from it the phase diagram for the system when we consider k to be a 
chemical potential rather than a Lagrange multiplier, that is, when we add — K$>[g]/T rather 
than —K,(<&[g] — Qq)/T to the action. Such a chemical potential for the Polyakov-Maldacena 
loop (PML) corresponds roughly to a density k of external quarks in the gauge theory. As in 
pure five-dimensional gravity analyzed at the end of the previous section, above the Hawking- 
Page temperature turning on k has no qualitative effect. Below that temperature, however, 
there is a second-order transition at k = to a phase with a tiny ten-dimensional hedgehog 
black hole. This black hole is smaller than the conventional ten-dimensional black hole, and 
exists only by virtue of the strings sourced by the external quarks, whose tension holds the 
horizon open. As k is increased, at some point there is a first-order transition in which the 
system flips to the large five-dimensional black hole. Unlike in the case of pure gravity (figure 
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||), the curve separating the tiny and large black hole phases cannot end, since the two phases 
have different symmetries (50(6) is broken in one and not the other). This phase diagram 
is sketched in figure ||[ Of course, this picture is only the simplest one compatible with what 
we know at present, and the true phase diagram could be more complicated. 

5.2 Gauge theory on R 3 



N 2 T 3 




Figure 9: Left: K versus <J> for (top to bottom) (2nT) 2 = 10, 20, 40, 80, oo. The bottom curve is the 
result for the gauge theory on R 3 . Note that each curve — including the one for R 3 — goes to a minimum 
value of <&. Here k is the string tension per unit volume on the boundary; /there — Right: 
The left figure blown up near the origin, showing the curves for (2ttT) 2 = oo (left) and 80 (right). The 
curve for T = oo presumably ends at the origin, representing AdS, although this region is difficult to 
access computationally. The solutions corresponding to the points a, b, c are plotted in figure [fl]. 

On R 3 , unlike on S s , different temperatures are related by conformal transformations, 
and therefore exhibit the same physics. There is therefore only one free energy diagram to 
compute, which is F/T 4 versus $/T. The only effect on the action ( |5.3| ) of going to R 3 is 
to remove the term 6e _2cr from the potential. The theory on R 3 can be considered as the 
T — > co limit of the theory on S" 3 . In figure ^ we plot k/T 3 versus $/T on 5 3 for a range of 
high temperatures, together with the result on R 3 . Three examples of hedgehog black brane 
solutions are plotted in figure 11. Figure 10 shows the free energy diagram on R 3 , obtained 



by integrating the (lowest) $— k curve in figure 
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Figure 10: Left: Free energy diagram for the gauge theory on R 3 . Here F is the free energy per 
unit volume on the boundary. Right: left figure blown up near the origin. 
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